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ABSTRACT 

To better constrain models of cool core galaxy cluster formation, we have used X-ray observations 
taken from the Chandra and ROSAT archives to examine the properties of cool core and non-cool 
core clusters, especially beyond the cluster cores. Using an optimized reduction process, we produced 
X-ray images, surface brightness profiles, and hardness ratio maps of 30 nearby rich Abell clusters 
(17 cool cores and 13 non-cool cores). We show that the use of double /3-modcls with cool core 
surface brightness profiles and single /3-modcls for non-cool core profiles yield statistically significant 
differences in the slopes (i.e., (3 values) of the outer surface brightness profiles, but similar cluster core 
radii, for the two types of clusters. Hardness ratio profiles as well as spectroscopically-fit temperatures 
suggest that non-cool core clusters are warmer than cool core clusters of comparable mass beyond the 
cluster cores. We compared the properties of these clusters with the results from analogously reduced 
simulations of 88 numerical clusters created by the AMR Enzo code. The simulated surface brightness 
profiles have steeper /3-model fits in the outer cluster regions for both cool cores and non-cool cores, 
suggesting additional ICM heating is required compared to observed cluster ICMs. Temperature and 
surface brightness profiles reveal that the simulated clusters arc over-cooled in their cores. As in 
the observations, however, simulated hardness ratio and temperature profiles indicate that non-cool 
core clusters are warmer than cool core clusters of comparable mass far beyond the cluster cores. 
T he general similari ties between observations and simulations support a model described in Paper 
I (|Burns et al.l l2008h suggesting that non-cool core clusters suffered early major mergers destroying 
nascent cool cores. Differences between simulations and observations will be used to motivate new 
approaches to feedback in subsequent numerical models. 

Subject headings: cooling flows, cosmology: observations, telescopes (Chandra, ROSAT), X-rays: 
galaxies: clusters 



1. INTRODUCTION 

Due to their high overdensity relative to the cosmic 
mean and their subsequent relative rarity, rich galaxy 
clusters have previously been treated as physically iso- 
lated objects (|Voitl l2005l . The combination of large mass 
and relative isolation make them an attractive tool for 
precisely measuring cosmological parameters s uch as r^ m , 
Qh and the dark energy equation o f state (e.g.lLin et all 
20031: ISanderson fc Ponmanl 120031: lLaRoaue et al.l 120061: 



Allen et al.l I200& IVikhlinin et al.ll2009T) . Assuming that 
the collapse of the cluster has captured all the mass 
within its gravitational sphere of influence, the relative 
densities of dark matter an d baryons should be the same 
as in the wider universe ([White et all 119931: lEke et all 
119981 : iFrenk et al.l ll999f). However, accurate determina- 
tion of the baryon fraction in the cluster, and thus for 
the universe as a whole, requires accurate measurements 
of both the baryon and dark matter mass of the system. 
Of the two types of galaxy clusters, cool core (hereafter 
CC) and non-cool core (hereafter NCC), CC clusters are 
more often chosen for these measurements because they 
are assumed to be dynamically relaxed and, thus, may 
be accurately fit with a density profile that provides the 
gas mass to dark matter ratio 



pronie tnat provid es tne 
(|Allen et all l200l . [20081: 
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IVikhlinin et"al|[2009l ). 

Cool core clusters demonstrate sharply peaked X-ray 
emission at their centers due to condensed regions 
of cooler gas that ar e brigh te r tha n the surrounding 
mater ial (see e.g. iFabianl (|1994f ): IDonahue fc Voitl 
(|200l for reviews of cool core clusters). In recent 
flux-limited sampl es, approximately half of the clusters 
contain cool cores (|Chen et al.ll2007h . While flux-limited 
samples are likely to increase the number of CC clusters 
found through the selection bias of a brighter core, CC 
clusters still clearly represent a significant fraction of 
galaxy clusters. T he lack of observed strong cooling 
flows in the se cores (jTamura et al.ll200ll : iPeterson et alJ 
120011 I2003D has been presented i n tandem with sev- 
eral recent AGN observations (lFabia"n et~aTl [20061 : 
iMcNamara fc Nulsenll2007l : IWise et al.ll2007h suggesting 
some form of AGN feedback inhibits run away cooling 
(iHeiirz et al.l 120061: iMathews et all l2006t iNusser et all 
20061 : ISiiacki fc Springell 120061 : Binnev et all 120071 : 



Cattaneo fc Tevssierll2007t ICiotti fc Ostri kcr 2007| ). The 



lack of runaway cooling can be interpreted as potential 
evidence that CC clusters have relaxed and achieved a 
relatively steady state of cooling and accretion making 
them more attractive for cosmological studies. 

One may find the baryon mass of galaxy clusters uti- 
lizing measured X-ray surface brightness radial profiles 
by fitting them with a /3-model, 

S x (r) = S x0 l+f , (1) 



2 
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where S x o is the normalization factor, r core is the size 
of the cluster 'core' and (3 is the slope of the power 
law. This parametric fi t of S x has been successful 
for some galaxy clus ters (|Jones fc Formanl [l98l fl999l : 
IVikhlinin et al.lTl999h . however, the fit implies no phys- 
ical processes or assumptions about the cluster. The (3- 
m odel was orig inally developed for an isothermal system 
bv lKind (|1966f ), and though clusters are known to be non- 
isothermal (e.g. iMarkevitch et all 119981 : IVikhlinin et alJ 
120051 12006D . the /3-model of ten yields reasonable empir- 
ical f its to observed profiles (jCavaliere fe Fusco-Femianol 
119761 ). Nevertheless, if a cluster is assumed to be isother- 
mal, then the X-ray surface brightness can directly trans- 
late into a function of gas density, 



n e (r) = n e0 1 + 



(2) 



Additionally, detailed Chandra and XMM observations 
show the dark matter potential i s consistent w i th the 
NFW profile para meterization of iNavarro et al.l (jl997f ) 
(jMushotzkvl I2004T ). The total cluster mass and rela- 
tive quantities of baryonic and dark matter constrain 
the values of i] m a nd ffy ([Wang fc SteinhardU 1 19981: 
Haiman et all 120011: lAllen et all l2008t IVikhlinin et all 



2009|) While/3-models have been shown to fit NCC clus- 



ters well, they are typically inadequate for CC clusters, 
owing to the exc ess surface brightness from their denser, 
cooler cores (e.g. IVikhlinin et al.ll2006h . To combat this, 
the over-dense CC cores are often excised from the j3- 
model fit, o r a second /3-model is included as well to 
fit the core dPratt fc Arnaudll2002t iPointecouteau et alJ 



l200l IVikhlinin et al.ll2006HSantos et al.ll2008h. 

In the light of r ecent work (e.g. iBurns et alJ 120081 : 
iJeltema et alJ 120081 ). two questions beg to be asked: 
What causes the X-ray properties of CC and NCC clus- 
ters to be so different? Are CC clusters as simple as 
previ ously assumed? The results of Paper I ([Burns et al.l 
120081 hereafter B08) predict some observed differences 
between CC and NCC clusters may result from different 
merger histories. Cosmological N-body/hydrodynamic 
adaptive mesh refinement simulations with cooling and 
feedback indicate that NCC clusters suffer major mergers 
earlier in their development than CC clusters, disrupting 
nascent cool cores. We defined a "major merger" in B08 
as a cluster encounter during which ~ 50% of the clus- 
ter mass is accreted over roughly 1 Gyr, a merger large 
enough to potentially destroy any existing cool cores. CC 
clusters, on the other hand, undergo only minor mergers 
so these cool cores are less susceptible to disruption. In 
the intermediate cluster region just beyond the cluster 
cool cores ( 0.05 — 0.3r2oo where r2oo is approximately 
the virial radius) the simulations find ~ 40% more cool 
gas in CC clusters than NCC clusters. This excess cool 
gas is also indicated by an offset between CC and NCC 
hardness ratio profiles. Additionally, single /3-modcl fits 
to surface brightness profiles in this intermediate region 
indicate that CC cluster flux is substantially overesti- 
mated there. This overestima tion leads to a signi ficant 
bias in estimated cluster mass (|Hallman et al.ll2006h , and 
suggests that CC clusters as cosmological probes must be 
used with caution. Motivated by these predictions, this 
paper looks for evidence of differing merger histories in 
observed clusters, further evidence that may cast doubt 



on the assumption that CC clusters are simple, relaxed 
objects. We also examine differences between observa- 
tions and our simulations in order to improve our clus- 
ter models. Using a sample of cluster observations from 
the Chandra and ROSAT archives, we look for evidence 
of differing merger histories between CC and NCC clus- 
ters in the hardness ratio profiles and spectroscopically 
fit temperatures of the observed clusters. We also com- 
pare observed surface brightness profiles and c alcula ted 
j3- model fits to the simulations of IBurns et al.l (|2008h . 

In §[2] we describe the observational cluster sample and 
the individual Chandra and ROSAT data used. We dis- 
cuss in § [3] the data reduction pipeline for both Chan- 
dra and ROSAT data, detailing differences between the 
standard pipeline and our process. We present obser- 
vational results including surface brightness profiles, /?- 
model fits, hardness ratio profiles, spectroscopically fit 
temperatures, and comparisons with published temper- 
ature maps from different X-ray observatories for vari- 
ous clusters in § 21 In § El we compare the simulations 
described in B08 with our observational results to deter- 
mine how our simulations may be improved in the future, 
and also find that the observed clusters indicate similar 
signs of different merger histories as found in our simula- 
tions. Finally, conclusions and a summary of the paper 
are in § [6l 

2. OBSERVATIONAL DATA SAMPLE 

Using both the Chandra and ROSAT archives, we 
chose 30 ri c h clus ters to analyze (17 CC and 13 NCC). 
lEdge et all (|1990f ) defines a sample of 55 clusters with 
fluxes greater than 1.7 xl0~ n ergs cm~ 2 s _1 as mea- 
sured with the Einstein or EXOSAT satellites, a reason- 
ably complete list of clust ers for these fluxes and above. 
We used a subset of the lEdge et aLI (jl990l ) list as our 
data sample; the sample consists of Abcll clusters with 
at least one observation of 3 ksec or greater from each of 
the Chandra, ROSAT and XMM archives. In the final 
sample, detailed in Tabled total Chandra exposures are 
all > 10 ksec and total ROSAT exposures are all > 3 
ksec. While the XMM data are not used in this paper, 
we plan to combine it with the ROSAT and Chandra 
results in a future publication. Clusters are included 
regardless of morphology or substructure present in an 
effort to limit selection bias. Given the diverse popula- 
tion of X-ray observations in the archives, our sample is 
not statistically complete, however we believe it is rep- 
resentative of CC and NCC clusters as described below. 
This is important if any significant conclusions about the 
properties and merger histories of general CC and NCC 
clusters are to be made. The clusters we analyzed are 
listed in Table [T] with their X-ray properties as well as 
Chandra and ROSAT observation details. 

Determinat i ons o f NCC or CC were taken from the 
lO'Hara et al.l (|2006h list that includes all of our clusters. 
However, we differe d from their list at one point. While 
lO'Hara et all (|2006t ) lists Abell 3562 as a CC cluster, we 
chang ed its designation to a NCC cluster. IChen et alJ 
(|2007h lists mass accretion into the core of the cluster as 
0.0 M Q /year and a cooling time of 1.3 x 10 10 years, both 
typical values for NCC clusters. The results we obtained 
from the data analysis also place it solidly into the NCC 
category. 

Rcdshifts were drawn from lO'Hara et all (|200l . while 
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M200 and r2oo were calculated from A/500 an d the tem- 
perature, respectively. Note that subscripts on M and r 
denote cluster mass and radius where the cluster overden- 
sities with respect to the critical density p cr equal these 
subscripts. Further note that the virial radius of the clus- 
ter is at ri78, which is ~ r2oo . For M2nn , we used the val- 
ues of M500 taken from lChen et all (120071) a nd the scaling 



relation M 200 = 1.377M Rnn fromlVoitl (12005T). For t he val 



ues of r 2 oo , we calculated r 50 o from the lChen et al.l (|2007h 
temp eratures using l the relation 7*500 = 0AA7h^n T2 527 



Mpc (jO'Hara et al. l2006t iFinoguenov et aLll2001[ ). The 



temperatures from lChen et al. M2007D were used in par- 
ticular as they are corrected for the biasing effects of 
cluster cool cores in average cluster temperature esti- 
mates. In reality, these Tchen temperatures are Th, the 
hot-co mponent cluster temperatures from llkebe et al.l 
(120021). We thus be lieve that the temperatures given 
bv IChen et all (120071) are more physically relevant esti- 
mates of overall cluster temperatures. Values of 7*500 were 
then scaled to 7*200 via r 2 oo = 1.51r5oo, which is based 
on an NFW dark matter density profile (jNavarro et alJ 
1997). The cosmological values assumed throughout this 
paper for the observations are Ho = 71 km s _1 Mpc -1 , 
fijvf = 0.27 and a "flat" universe where f2 = l 3 . 

Figures [1] and [5] demonstrate the range of clusters in 
our sample and hopefully further demonstrate that our 
sample is representative of CC and NCC clusters, if not 
a statistically complete set. If clusters are scale-free, one 
should find a correlation between temperature and mass 
(Kaiser 1986), which is clearly visible for our sample in 
Figure [TJ An empirical mass-temperature relation for 
a large sample of CC and NCC clusters was found by 
IChen et all pool to be 



logi 



Mi 



.Mill 



5 x 10 14 M 







-0.112+ 1.54 logj 



T 



4keV 



(3) 

This relation is over-plotted in Figure [TJ and shows that 
our clusters contain quite representative masses and tem- 
peratures from the broader sample. 

As noted in B08 and seen in the data from 

10'Hara et all ll200l and IChen et ail (|200l . CC clusters 
are somewhat segregated from NCC clusters in tempera- 
ture and mass as we see in Figures [TJ and __l In both ob- 
served and simulated samples, NCC clusters tend to have 
higher overall temperatures and mass. The CC clusters, 
on the other hand, are observed to fill the whole range 
of temperatures and mass in the sample but have some- 
what lower mass and are cooler on average than NCC 
clusters. The fact that NCC clusters are only found with 
relatively high mass and temperatures corroborates our 
suggestions that CC and NCC clusters experience differ- 
ent merger histories as discussed in B08. The segregation 
in cluster types is fully explained if NCC clusters expe- 
rience major mergers during their formation, inhibiting 
their ability to retain cool cores and warming their ICM 
in the merger process, while CC clusters grow via only 
relatively minor mergers later in their history once their 
cool cores have been established. 



3 Note that the cosmological values used in Chen et al. (2007) 
are Ho = 50km s — 1 Mpc -1 , Q^f = 1, and = 0- The masses 
quoted in Table [T] and used throughout this work are corrected to 
our cosmology. 



With the exception of Abcll 2204, the sample consists 
of nearby rich clusters with redshifts between 0.01 < z < 
0.09. In Figure [H one notices that the (related) cluster 
distributions in temperature and mass versus rcdshift are 
reasonably random. A slight correlation, though, is vis- 
ibly present in that the NCC clusters are clumped at 
higher rcdshift and temperature. This correlation is not 
necessarily a selection bias, but an expected consequence 
of the cosmological volume element. At higher redshifts, 
the larger cosmological volume and larger numbers of 
clusters yields more of the rarer high mass clusters (NCC 
clusters). 

Though statistically incomplete, our Chandra / ROS AT 
X-ray subsample appears broadly representative of the 
X-ray pr operties of the lar ger statistical sample of clus- 
ters from lEdge et al.l (|1990f ) and expectations from mod- 
els. The temperature-mass correlation in our subsample 
follows the empiri cal relation f or ge neral CC and NCC 
clusters found by IChen et all (|2007l ) . The clusters are 
also distributed as expected in rcdshift space. Finally, 
the slight correlation of mass and temperature with red- 
shift observed in our subsample is in accordance with 
effects from the cosmological volume element. We feel 
our subsample is representative enough of a larger, more 
statistically complete cluster set to draw general conclu- 
sions from comparisons with simulated clusters in B08. 




T (keV) 

Fig. 1. — Plot of cluster mass versus temperature for the sample 
in Table[T] Red points correspond to NCC clusters and blue points 
correspond to CC clusters. The mass-temperature relation from 
IChen et al.l (t2007h is over-plotted, showing that our cluster sample 
is similar to a more statistically complete cluster set. 



3. X-RAY DATA ANALYSIS 

3.1. Chandra Data Reduction 

While reducing the Chandra data, we carefully fol- 
lowed the well-tested procedures outlined in the CIAO 
ACIS data pipeline threads 4 and the Markevitch cook- 
book for ACIS background datasets 5 . In order to ensure 
that all of the data used were reduced accurately with 
the most recent calibration files in a consistent man- 
ner, we started observational reductions with the raw 
Level 1 telemetry files and reset all status bits associated 
with each event. Using CIAO (version 3.3.0), we calcu- 
lated new "bad pixel" maps with the acis_run_hotpix 
tool for each individual observation. Unlike the stan- 
dard data pipeline, however, we did not discard the 

4 http:/ /cxc. harvard.edu/ciao/guides/acis_data. html 

5 http:/ /cxc. harvard.edu/contrib/maxim/acisbg/COOKBOOK 
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Fig. 2. — Plots of cluster temperature and mass versus redshift 
for the sample in Table [T] Red points correspond to NCC clusters 
and blue points correspond to CC clusters. A weak correlation in 
mass and temperature with redshift is detectable, but this is an 
expected consequence of the cosmological volume element. 

colu mns next to bad p i xels, consistent with recent ad- 
vice l|Markevitch et all I2003D . Using the CIAO tool 
acis_process_events, we corrected each observation for 
differing gains across the CCDs, time-dependent gain, 
and CTIs (charge transfer inefficiencies) to recreate a 
Level 1 file using the latest calibration tools and files. 
When an observation included chip S4, wc utilized 
the tool destreak to remove systematically hot pixel 
columns that occasionally appear in tha t chip. Events 
with unacceptable ASCA-based grades (jBaganoff et alJ 
I2003T) and point sources were discarded. We also used 
the acis_process_events tool to apply a randomization 
to both the position of the event within the given pixel 
that detected the photon as well as the recorded energy 
within the PHA (Pulse Height Amplitude) bins. This 
randomization removes stratification of the data at fine 
position and energy bins. 

To remove unresolved flares in the same manner as is 
used for the ACIS background datasets, event data for 
front- and back-illuminated chips were first separated. 
Good Time Intervals (GTI) were determined using a non- 
CIAO tool called lc_clean, aga in recommended and writ- 
ten bv lMarkevitch et ail (|2003t ). This resulted in excised 
periods of time where the incident flux differed from the 
observational average by more than 20%. The cleaned 
event files were then merged, resulting in the new Level 
2 event files for each cluster observation. 

3.1.1. Background Subtraction 

Due to the extended spatial structure of the clusters 
and the limited field of view of Chandra, using an exte- 
rior portion of each image as a representative background 
is inadequate. Instead, we used the ACIS background 
files discussed in the Markevitch cookbook. We followed 
the procedures laid out in the CIAO threads 6 in order 
to appropriately match the backgrounds to each clus- 

6 http:/ /cxc. harvard.edu/ciao/threads/acisbackground/ 



ter observation. Most of our cluster observations span 
multiple chips, so the background files, which are sepa- 
rated by chip, were then merged while maintaining sepa- 
rate files for front- and back-illuminated chip sets. If the 
processing of the background and event files used differ- 
ent gain files, then the backgrounds were re-calibrated 
with the CIAO tool acis_process_events to match the 
gain files used on the observation event files. Ideally, 
the calibration files used for CTIs and time-dependent 
gains would also match between our event files and the 
ACIS background files, though it is not possible to up- 
date these gains. According to the CIAO threads, how- 
ever, the induced error is expected to be minimal for 
extended sources like the clusters in our observations. 
Lastly, we re-projected the merged backgrounds into the 
coordinates of the observations using the tool repro- 
ject_events and the aspect solution files from the ob- 
servations. 

To determine the exposure duration of the backgrounds 
in order to normalize them to each observation, the im- 
age GTI duration was multiplied by the ratio of the 9- 
12 keV flux between the image and background files. In 
this band Chandra's effective area is essentially zero, thus 
nearly all the counts are particle background. The obser- 
vations and backgrounds were then split into two energy 
bands (0.5-2.0 keV and 2.0-8.0 keV) and each band was 
thereafter treated as a separate observation in order to 
construct hardness ratio maps later in the analysis. Mul- 
tiple observations of the same cluster were co-added into 
a single image and backgrounds were subtracted with the 
tool dmimgcalc. These background-subtracted mosaic 
images constitute the final "counts" images and were used 
for the hardness ratio analysis described below. 

3.1.2. Exposure Maps 

Each Chandra observation was taken over an extended 
period of time (tens of kiloseconds) and the pointing dur- 
ing that time both drifts and is dithered to prevent pixel 
differences from biasing the image. While pointing is 
corrected within the Level 1 file, aspect files from the 
pointing history were calculated to help determine the 
effective area convolved across the CCD surface. For 
this correction, we combined the spatial region of the 
cluster (minus the cool core (< 0.05r2oo) when appro- 
priate) and extracted the observational spectrum. This 
was t hen fit with a combination of g alactic Nh absorp- 
tion orrison fc McCammonl I1983D and a Ra ymond- 
Smith plasma model ( Raymond fc Smith! Il977f ) within 
the Sherpa program to recreate the expected spectrum 
reaching the front of the telescope. This spectrum, the 
instrument responses (CIAO's ARF and RMF files) and 
observational GTIs were then folded together as detailed 
in the standard CIAO pipeline to create an exposure map 
for each observation, containing effective area, quantum 
efficiency and pointing. We then divided the exposure 
maps from the counts images to give a final product in 
units of surface brightness (photons s _1 cm _2 arcsec -2 ). 

These exposure maps were used only for surface bright- 
ness profiles (0.5-2 keV) and not in the hardness ratio im- 
ages (0.5-2 keV and 2-8 keV) as we discuss in Section 2J 
Because of the wide energy bands and spatial coverage, 
there is no exposure map that can accurately represent 
the effective area convolved with the source spectrum, es- 
pecially over the 2-8 keV range. The lower energy range 
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of 0.5-2 keV has a relatively flat effective area function 
and the spatial change across the field of view utilized is 
on the order of 1%. Thus, while the exposure maps were 
used to calculate fluxes in this band for surface bright- 
nesses, we followed a different approach for the images 
used in the hardness ratios. 

Figure [3] shows four examples of reduced Chandra flux 
images, two CC and two NCC clusters. The final mo- 
saicked images have square-root scaling and cover the 
soft band 0.5-2.0 keV. The images were also smoothed 
with a three pixel FWHM gaussian. The flux images 
were binned by a factor of eight so that one pixel corre- 
sponds to 12 arcsec. The CC clusters are very distinct 
from the NCC clusters, with significantly lower flux in 
the outer cluster extent in contrast with a sharp flux rise 
in their cores. The NCC surface brightnesses exhibit a 
much shallower drop-off moving outward from the clus- 
ter centers. The images for A478, A1795, and A3667 are 
mosaics of three individual pointings, while the image of 
A401 is made from two. 

3.1.3. CIAO Calibration Updates 

Since the original work in this analysis was completed 
with CIAO 3.3.0, a newer version of CIAO and the 
CALDB (both versions 4.1.1) have been released. The 
new calibration files take into account a contaminant 
originally missed during ground calibration tests that 
helps explain both discrepancies in high-temperature 
clusters between Chandra and XMM-Newton observa- 
tions and between the Chandra ACIS-I and ACIS-S in- 
struments themselves 7 . To make sure the use of the new 
calibration files would not significantly alter our findings, 
we re-analyzed six clusters with temperature > 6 keV 
(A85, A2255, A478, A401, A2319, A3266). We found 
changes of less than 5% in derived hardness ratios (see 
Section 14. 2p for all clusters other than A478, which it- 
self saw only 5 — 10% deviations. Derived temperatures 
for each of the clusters in our temperature analysis be- 
low (A85, A2255, A478, and A401; see Section _OJ| all 
dropped slightly compared to the results using CIAO 
3.3.0. Because the changes in hardness ratios were so 
small, the results discussed below in Section 14.21 derive 
from the analysis with CIAO 3.3.0 and its corresponding 
CALDB. The findings in Section [43] for clusters A85, 
A2255, A478, and A401, however, arise from the use of 
CIAO and CALDB 4.1.1 to take advantage of the new 
calibrations. 

3.2. ROSAT Data Reduction 

While Chandra is the superior instrument in terms of 
effective area, calibrations, and breadth of energy re- 
sponse, the ROSAT PSPC has a field of view approx- 
imately four times larger. It is important to note, how- 
ever, that ROSAT observations in our cluster subsam- 
ple generally have much shorter exposure times than 
the Chandra observations of the same cluster. Data 
we re reduced usin g the process and tools as described 
bv ISnowden et all (|1994f ). and all the ROSAT pointed 
observations in our subsample used the PSPC detector. 
ROSAT data are not available at an uncorrected level 

7 http:/ /cxc. harvard.edu/caldb/downloads/Release_notes/ 

CALDB_v4.Ll.html 



equivalent to Chandra data. Thus, the event files pro- 
vided were fully corrected for most inconsistencies and 
we did not recalibrate the events other than the screen- 
ing discussed below. Using the same tools and limits 
as the Chandra pipeline, the time delineated flux was 
screened for flare events (time periods with fluxes more 
than 20% from the observational mean), which were re- 
moved from the observation. Events also determined to 
be "spurious trailing" events (recorded events that repre- 
sent the residual charge from an actual prior photon) and 
point sources were excised. The events were binned by a 
factor of 30 to form 15 arcsec x 15 arcsec pixels with the 
goal of increasing the S/N, and finally cut to an energy 
range of 0.42-2.01 keV. 

3.2.1. Exposure Maps 

Unlike Chandra, the ROSAT data reduction pipeline 
does not consider the spectrum of the source when cre- 
ating an exposure map. While this method is inher- 
ently inaccurate, the changes in both the spectrum and 
ROSAT's effective area across the 0.42-2.01 keV energy 
band are small enough to make the result acceptable. 
The pointing of the telescope during the exposure and 
the GTI time were combined with the effective area and 
quantum efficiency calibrations. These effects are con- 
tained in a previously calibrated "detector map" provided 
by the Snowden package. Since the closest energy band 
of the detector maps is 0.42-2.01 keV, the ROSAT and 
Chandra images do not exactly match in energy Counts 
images were finally combined into a mosaic and corrected 
by a similar mosaic of exposure maps to create a final im- 
age. 

3.2.2. Background 

Unlike the Chandra images, ROSAT has a sufficient 
field of view to use the outer regions of an image as a 
background source. A flux average of the annulus at the 
exterior of an image was subtracted from the rest of the 
image to create the final background-subtracted prod- 
uct. This annulus varied in radius between images due 
to detector edge blurring and noise. This was typically 
20 pixels in width and ended a few pixels short of the 
largest radii deemed free of these edge effects. 

4. OBSERVATIONAL RESULTS 

Numerical cluster simulations from B08 suggest ob- 
servable differences in the X-ray properties between CC 
and NCC clusters, owing to different evolutionary path- 
ways. In particular, our numerical models predict differ- 
ences in surface brightness and temperature/hardness ra- 
tio profiles in the region beyond the cluster cores. To test 
these predictions, we studied surface brightness, hard- 
ness ratio, and temperature profiles derived from our 
reduced Chandra and ROSAT observations, which are 
described in detail below. In SjS] we compare analogous 
simulated profiles to these observed results to draw con- 
clusions not only about the merger histories of observed 
CC and NCC clusters, but also about what physics and 
feedback parameters need further development in next- 
generation simulations. 

4.1. Radial Surface Brightness Profiles 

We created flux images (photons s _1 cm _2 arcsec -2 ) for 
both Chandra (0.5-2.0 keV) and ROSAT (0.42-2.01 keV) 
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Fig. 3. — Four example Chandra flux images of CC and NCC 
clusters over the soft band of 0.5-2.0 keV. Each is the result of 
multiple individual observations that have been combined to form 
a final data product and then smoothed with a 3 pixel (~ 12 arcsec) 
FWHM gaussian. The FWHM for A478, A1795, A401, and A3667, 
respectively, are the following: 19.6 kpc, 14.2 kpc, 16.9 kpc, 12.2 
kpc. Each image is cut to a radial metric field-of-view of 0.3r2oo- 
All four images are square-root scaled and show — 2 X 10 — 5 (CC) 
and — 2 X 10 — 6 (NCC) photons s _1 cm _2 arcsec — 2 , respectively. 

observations from the reduced counts images by perform- 
ing exposure map corrections, as detailed above. The 
point sources, as well as obvious instrument artifacts such 
as CCD boundaries, were masked out in each image be- 
fore taking averages of radial annuli. While real galaxy 
clusters are not perfectly symmetric azimuthally, the dif- 
ferences are usually minimal even for clusters that have 
obvious structure and an annular average is a reasonable 
repr esentation of the stat e of the cluster at that radius 
(e.g. IVikhlinin et al.|[l999f) . The size of radial bins were 
adaptively binned for each cluster separately to maxi- 
mize both the number of bins and S/N. Figure [4] shows 
the combined Chandra and ROSAT surface brightness 
profiles for all clusters in our subsample. 

Several traits are common across most of the observa- 
tions as shown in Figured] First, the ROSAT data were 
included in this work because it was hoped that with its 
larger field of view the surface brightness profiles could 
be extended to larger radii. This was largely untrue. The 
errors of the ROSAT data are often large enough to pre- 
clude a confident fit even where the Chandra data have 
low errors. While some nearby clusters do gain from the 
ROSAT inclusion, on the whole it docs not add signif- 
icantly to the results. As a consequence, ROSAT data 
were excluded from cluster profiles for which Chandra 
data extended to roughly 0.4r2oo- Second, the reliable 
ROSAT data and the Chandra data match well across 
every cluster but two, A426 and A2063, so ROSAT data 
are excluded for these clusters as well. It is unknown why 
the data for A2063 do not agree as well as in the other 
clusters. Chandra data for A426 show obvious substruc- 
ture that does not appear in the ROSAT data, leading to 
deviations in the radial profile averages that account for 
the poorer match. In the end, ROSAT data were used for 
A85, A119, A262, A1060, A1367, A1656, A2199, A2319, 



and A4038. Finally, one may notice that 0.3 — 0.5 r2oo is 
a reasonable median of the outer radius at which clusters 
can be reliably fit with a /3-model based on the standard 
deviations of values at each radial bin. NCC clusters 
are fit with a single /3-mod el, while CC cluste rs are fit 
with a double /3-model (e.g. iSantos et al.ll2008ft . The /?- 
modcl fits are done using the mpfit routine, which uses 
a Levenbe rg-Marquardt leas t -squares alg orithm for min- 
imization (iMarkwardtl 12001 IMordli977l) . The use of a 
double /3-model is to account for the influence of the cool 
cores, which may contaminate the surface brightness pro- 
files beyond the cores themselves, potentially biasing a 
single /3-model fit that only excludes the fiducial cores. 

In Figure [51 the averaged surface brightness profiles of 
30 clusters separated by cluster type are shown with the 
average of /3-model fits (dashed lines). The averaging was 
done to increase S/N in the outer profiles. To create the 

averaged profiles, each cluster was first scaled by M^q 
to account for l uminosity differences due to the differing 
cluster masses (|Finoguenov et all 120011 ). The rescalcd 
profiles were radially binned between fixed fractions of 
r2oo to standardize the data, and these binned cluster 
fluxes were then averaged to produce the data in Figure 
Uncertainties are shown as 90% regions, demarcated 
by dotted lines, within which 90% of our data can be 
found. Each cluster was fit individually with a /3-modcl 
between and 0.4r 2 oo- The results of the individual clus- 
ter fits were then averaged to obtain an average r core and 
(3 values in the case of NCC clusters or two r core and two 
/3 values in the case of CC clusters. The resulting av- 
erage /3-model fits are overplotted in Figure [5J The size 
of deviations from the mean increase past ~ 0.2 t^oo in 
CC clusters, demonstrating that it is difficult to assess 
the validity of a /3-modcl fit to the outer region of these 
clusters. As is evident in Figure|31 the surface brightness 
profiles for many CC clusters do not go out as far as those 
for NCC clusters, which is contributing to the larger un- 
certainties in the outer CC profile. This is likely to be 
a selection effect as NCC clusters are in general farther 
away and so cover a smaller angular area, which results 
in better sampled surface brightness profiles farther from 
the center of the cluster. However, this does not exclude 
other potential reasons for larger deviations in the outer 
average CC profile. Perhaps some physical effect has a 
role. Are the CC clusters truncated compared to NCC 
clusters? Are NCC clusters being "puffed out" by massive 
mergers and heating? Further investigation is needed to 
explain where the larger deviations from the mean outer 
CC surface brightness profiles may be originating. 

Some differences in surface brightness profiles between 
CC and NCC clusters are evident in Figures [4] and [5j CC 
clusters have considerably higher emission in their cores 
due to the cooler, denser gas. NCC clusters have much 
flatter emission in the inner cluster region that begins to 
drop off at ~ 0.1r2oo- The overall shapes/slopes of the 
profiles beyond the cores appear to be different for these 
two cluster types. To help quantify these differences, 
Figure [6] shows histograms of the /3-modcl parameter fits 
for each cluster sample. The averages of the individual 
parameter fits are given in Table [21 and the individual 
best-fit values and fit statistics are presented in Table 
[31 T wo-sample, one-tailed Kolmogorov-Smirnov (K-S) 
tests ([Wall fc Jenkinsl[2003T) were also performed to mea- 
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sure the probabilities that CC and NCC model fits come 
from the same parent distributions. We draw two impor- 
tant conclusions from the data in Figure [H and Table [21 
The first is that the outer slopes of the surface brightness 
profiles are statistically different between CC and NCC 
clusters, with a probability they come from the same dis- 
tribution of only 0.002. Could this be caused by the dif- 
ferent evolution history of CC and NCC clusters? This 
will be discussed further in The second conclusion 
we draw is that the r core values are statistically indistin- 
guishable {Prob = 0.209) between CC and NCC clusters 
This differs from previous results (jVikhlinin et al.lf2006l : 
iBurns et al.ll2008f ) where single /3-models are used for CC 
clusters. One expects the average r core value for CC and 
NCC clusters to be similar if the data were completely 
corrected for the cooler cores. If cool cores are not com- 
pletely removed, however, then single /3-model fits to CC 
clusters may bias the r core results low, resulting in dif- 
ferent r core results for the two cluster types. 



A262 (CC) 



A426 (CC) 



</<m 
A401 (NCC) 



A478 (CC) 



Fig. 4. — Observational surface brightness profiles with both 
Chandra (red) and ROSAT (black) data. Error bars are the stan- 
dard deviation of each data bin from their respective telescope. 
The ROSAT error bars are represented by the black dotted lines 
for clarity purposes and disappear when the measurement minus 
the error becomes negative. Black dashed lines show the /3-model 
fits from to 0.4 r^oo f° r each cluster. A double /3-model was used 
for CC clusters to account for the presence of the cores, while sin- 
gle /3-model fits were used for NCC clusters. Vertical dashed lines 
at 0.4 r200 show the extent of data used for the /3-model fits. The 
double /3-model fits for CC clusters agree well with the observed 
data, except for A2204 where t here is reason to believ e the Chandra 
calibration data is suspect (sec Rciprich et al. 2008). 



4.2. Hardness Ratio Maps 

Motivated by the predictions of simulations from B08, 
we looked for evidence of differences in cluster evolution 
between CC and NCC clusters by exploring the region 
beyond cool cores to ~ 0.3 r2oo- The simulations indicate 
that CC clusters contain ~ 40% more cool gas in the re- 
gion beyond ~ O.f r2oo than NCC clusters. One way to 
estimate the cluster temperatures in this region is with 



the use of hardness ratio maps. A two-band hardness 
ratio map allows one to probe farther from the cluster 
core since it requires fewer counts than spectral temper- 
atures for a given signal-to-noise. The hardness ratio of a 
cluster's X-ray emission is also a proxy for temperature 
because a significant component of cluster X-ray emis- 
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sion is thermal bremsstrahlung radiation. Recall that 
for thermal bremsstrahlung 

S x <x [ n 2 e K{E,T)dE oc f ' nlT- 1/2 C - E/kT G(E,T)dE, 

J J E x 

(4) 

where S x is the X-ray luminosity in the energy band 
E\ — E2, G(E,T) is the quantum correction G aunt fac- 
tor, a nd n e is the electron number density (e.g iPeacockl 
I1999T ) . From the temperature dependence of S x , one sees 
that the hardness ratio S X 2/ ' S x \ is just related to the ra- 
tio of temperatures T2 and T\ , which after integration is 
roughly proportional to (T2/T1) 1 / 2 . Thus, a higher hard- 
ness ratio corresponds to harder X-ray emission, which 
in turn implies a higher temperature. 

The X-ray spectrum of a cluster, however, has other 
contributing factors besides thermal bremsstrahlung con- 
tinuum emission. Line emission is an important spectral 
component, especially if the cluster is relatively cool. To 
demonstrate that the hardness ratio is in fact propor- 
tional to cluster temperature, we used the X-ray spec- 
tral fitting package XSPEC to model spectra of clusters 
of increasing temperatures from 1 to 10 keV. Metalli- 
cally was assumed to be 0.3 Z Q . The modeled spectra 
were broken into two standard Chandra bands, 0.5 — 2.0 
keV and 2.0 — 8.0 keV, the photon flux in these bands 
was integrated, and the ratio of hard band to soft band 
was calculated. The resulting temperature dependence 
of the hardness ratios of the modeled spectra is plotted 
in Figure [7] The temperature itself has little dependence 
on metallicity; we found only a few percent change in 
hardness ratio calculated in this way when varying the 
abundance between 0.3 and 0.7 Z Q . The square root de- 
pendence of the hardness ratio on the temperature for 
a bremsstrahlung model is also plotted. It is clear that 

A3266 (NCC) A3391 (K'CC) 




10 6 
10 5 

x 10" 







- 


Cool Core 




N^n^-t^o^^o^ 


_ Data 




Fit 





0.01 0.10 

r / r 200 

Fig. 5. — Average surface brightness of 15 CC (blue) and 13 
NCC (red) clusters observed with Chandra and ROSAT. (The j3- 
model fits for A426 and A478 are far outliers and are not included 
here.) Each cluster is scaled by Mjqq to account for luminosity 
differences due to mass. Dotted lines delineate the 90% regions 
of each data bin. The NCC clusters are arbitrarily divided by 10 
to separate the curves for clarity. Dashed lines are the averaged 
individual beta model fits out to 0.4r200- A double /3-model was 
used for CC clusters, while a single /3-model was used for NCC 
clusters. 
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Fig. 6. — Histograms of /3-model fit parameter values. Red is 
for NCC clusters and blue for CC clusters. (Top) Outer r core fits. 
(Bottom) {} fits. The fits for A426 and A478 are far outliers, and 
are not included here. It is clear from the histograms that CC 
and NCC r core fits are nearly identical, while there is a significant 
difference in /3 slopes. 

hardness ratio does increase with temperature, but not 
as fast as pure bremsstrahlung emission predicts. 

The hardness ratio images for the observed cluster sub- 
sample were created by dividing two Chandra images, the 
2.0-8.0 keV bandpass image by the 0.5-2.0 keV bandpass 
image. These images are counts-based and reduced as 
described previously. They lack corrections for the ef- 
fective area, quantum efficiency, pointing deviations and 
exposure durations, which are corrected by other means 
described below. 

First, the actual exposure time is inconsequential as 
each bandpass image was drawn from the same observa- 
tion and share the same exposure duration. As pointing 
was corrected in each event when recorded, the exposure 
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FIG. 7.— Hardness ratio (2.0-8.0 keV / 0.5-2.0 keV) for 
modeled clusters of temperatures 1 — 10 keV, including thermal 
bremsstrahlung and line emission. Also plotted is the expected 
hardness ratio including only thermal bremsstrahlung. The hard- 
ness ratio clearly increases with temperature, but not as quickly as 
expected for pure bremsstrahlung emission. 

map accounts for the changes of effective area as the tele- 
scope was dithered across the sky. This means that in 
theory each pixel could have a different effective area. 
We took the hardness ratio maps, however, out to only 
~ 0.3r2oo (~ 500 kpc), a range over which the effective 
area of the detector is reasonably flat and where there 
are sufficient counts when the cluster is centered in the 
field of view. 

We next accounted for the absolute difference between 
the energy bands. The difference in quantum efficiency 
and effective area between the two energy bands was re- 
moved by normalizing each to the idealized spectrum of 
the cluster. This idealized spectrum was found by fold- 
ing the instrument response into the expected spectrum 
of the cluster at the temperature given by Equation 03 
an Nh absorption model, and an APEC plasma model 
with 0.3 Z Q . Using WebSpcc 8 , an internet interface ver- 
sion of XSPEC, we combined the instrument response 
with the APEC model spectrum, providing the expected 
count rate in each energy band for gas at the M — T 
average temperature from Equation [3J The ratio of the 
two fluxes gives the hardness ratio of gas at this average 
temperature; each cluster's hardness ratio map was then 
normalized to this ratio. Thus, if we assume the aver- 
age temperature is in fact the virial temperature for the 
cluster, a hardness ratio of unity would correspond to a 
pixel at the virial temperature; less than one would be 
cooler, and greater than one would be warmer. 

Before division, we adaptively smoothed the images 
with a gaussian kernel using the CIAO tool csmooth. 
The gaussian size scales at each pixel using a Poisson 
expectation of noise to create an image with a signal-to- 
noise ratio of 3 in each pixel within the 0.5-2.0 keV band. 
The same smoothing kernel was applied to the high en- 
ergy band, the two images were divided, and the resulting 
hardness ratio image was normalized to the average clus- 
ter temperature using the WebSpec expected flux ratio. 
The final image gives an approximate visual representa- 
tion of the temperature across the cluster with respect 
to the virial temperature. This is directly comparable 
to the hardness ratio maps created from the simulations, 
which are also normalized to their virial temperatures. 

8 http:/ /heasarc. gsfc.nasa.gov/webspec/webspec. html 



Figure [8] shows examples of hardness ratio maps for six 
observed clusters. 

As a check on this process and in order to probe the 
histories of CC and NCC clusters, we examined clusters 
paired as close as possible in mass. This also translates 
into similar virial temperatures according to the mass- 
temperature relation (Figure [IJ . In this way, observed 
differences in the clusters cannot be attributed to differ- 
ent cluster masses. Matching mass to within 10%, we 
found three pairs of observed CC and NCC clusters from 
the subset of our sample for which hardness ratio maps 
were produced (see below for more details). One pair 
included A2204 and so was disregarded. The remaining 
pairs, along with one more unmatched pair for illustra- 
tion, and their bulk temperature fits from £ )4.3I for the 
region in question are found in Table [6j 

Not every cluster has sufficient Chandra coverage to 
calculate hardness ratios out to our standard 0.3r2oo- 
As a result, we used only 8 CC and fO NCC clusters in 
the hardness ratio data analysis. We also excluded Abell 
3571 and 2597 from the analysis even tho ugh they have 
suffici ent coverage. While both we and lO'Hara et al.l 
(|200l label A3571 as a CC cluster, iChen et al.1 (Ejol f 



lists its cooling time and mass inflow rate as 0.84 x 10 yr 
and 35 M Q / yr respectively. The cooling time is high and 
the mass inflow rate is low compared to other CC clus- 
ters. It also falls as an extreme outlier compared to the 
other CC clusters' HR values in our analysis. It is our 
belief that A3571 is best described as an "in-between" 
cluster, not an NCC, but not fully a CC either. Thus, 
we disregarded it for the HR analysis so as not to signif- 
icantly skew the CC HR mean. A2597 is also a peculiar 
CC outlier, with a rapidly increasing hardness ratio pro- 
file, which matches neither the typical CC nor the typical 
NCC result. Again, to avoid confusion from potentially 
including an "in-between" cluster, A2597 was removed 
from the hardness ratio analysis. 

The top of Figure [5] shows the individual normalized 
hardness ratio profiles for the observed CC (blue) and 
NCC (red) clusters. Standard deviation error bars for 
A3391, which are representative of typical errors in all 
the cluster profiles, are included for illustration. With 
the exception of the innermost regions of the CC clus- 
ters, the hardness ratio profiles for all clusters tend to 
be relatively flat out to ~ 0.3r2oo despite the modest 
decline in temperature with radius (see £|4.3p . This flat- 
ness can be attributable to the hardness ratio's modest 
dependence on temperature. In general, NCC clusters 
have higher normalized hardness ratio profiles than the 
CC clusters. To quantify this offset, we averaged the 
CC and NCC profiles, and these averaged profiles can be 
found in the bottom of Figure O with standard deviation 
error bars. Though the scatter is somewhat large, there 
is an average offset in hardness ratio well beyond the cool 
cores between CC and NCC clusters, with NCC clusters 
having larger ratios by ~ 0.1. 

To determine if the observed hardness ratio offsets 
are statistically signifi cant, we performed a K-S test 
(|Wall fc Jenkinsl l2003t ) on the unbinned hardness ratio 
map pixels from all cluster hardness ratio maps. We 
separated the sample sets into CC and NCC groups, 
and constricted the hardness ratio maps to the region 
(0.15 — 0.3)r2oo to be certain we did not bias the distri- 
butions by including any cool core pixels. Cumulative 
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distribution functions of the individual hardness ratio 
map pixel values are displayed in Figure \W\ where it 
is clear there is a significant offset between CC and NCC 
hardness ratios. Including unbinncd map pixels, the two 
sample sizes in the K-S test are quite large. As a result, 
the significance of the K-S D-statistic follows the signif- 
icance of a x 2 distrib ution with two degrees of freedom 
l|Wall fc Jenkinsl2003| ). The x 2 statistic corresponding to 
the obtained D-statistic is \ 2 = 29163, which means the 
probability that the CC and NCC hardness ratio pixel 
sets come from the same distribution is <gC 1%. We con- 
clude that the hardness ratio offset seen in the averaged 
profiles results from intrinsic differences in the properties 
of clusters of very similar mass. Given that hardness ra- 
tios track overall cluster temperature, this suggests that 
CC clusters have more cool gas than NCC clusters well 
beyond their cool cores. 




Fig. 8. — Six observational hardness ratio maps. Each is the two 
Chandra images, 2-8 keV/0.5-2 keV, normalized by the hardness 
ratio corresponding to the cluster virial temperature. White circles 
are locations of removed point sources. Each image is cut to a 
radial metric field-of-view of 0.3r200- (Top) CC clusters A478, 
A2204, A2029. (Bottom) NCC clusters A2255, A3158, and A3391. 
Note the clear difference in normalized hardness ratio between CC 
and NCC clusters. 



4.3. Bulk Temperatures 

The results of the hardness ratio analysis suggest 
that the offset seen in hardness ratios should also be 
present in mass-matched cluster temperature profiles. To 
test this hypothesis, we defined four regions for a sub- 
sample of mass-matched clusters evenly spaced between 
~ (0.1 — 0.3)r2oo and found bulk temperatures for each of 
these regions by fitting their X-ray spectra using XSPEC 
v. 12.3.0. The cores of CC clusters are much cooler than 
the outer cluster regions, and similarly, the central gas 
in an NCC cluster is often warmer than gas farther out. 
To ensure we removed all the central gas so as not to 
bias our cluster bulk temperatures low or high for CC 
or NCC clusters, respectively, we looked at the surface 
brightness profiles ( fc|4.1j) of the CC clusters in question, 
found where they begin to turn over at the edges of the 
cool cores, and doubled this radius, which was then used 
as an inner radius cutoff. These cutoffs were indepen- 
dently defined for each CC cluster, but have an average 
of 0.11 r2oo, which we used as the inner radius cutoff for 
all NCC clusters. 

We removed point sources with the CIAO tool wavde- 
tect, searching over five scales, 1, 2, 4, 8, and 16. With 
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Fig. 9. — (Top) Individual observed normalized hardness ratio 
radial profiles. Blue lines are CC clusters, whereas red lines are 
NCC clusters. Standard deviation errors are plotted for one cluster 
that are representative of all the clusters. (Bottom) Averaged CC 
and NCC hardness ratio profiles, with standard deviation error 
bars. Here and throughout the paper, CC and NCC clusters have 
been slightly offset radially for clarity. 
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Fig. 10. — Cumulative distribution functions for all independent 
and unbinned observed CC (blue) and NCC (red) hardness ratio 
map pixels in the range 0.15 — 0.3 T200 ■ While the distributions are 
quantifiably very different, as evident in £14.21 the distributions are 
quite visibly different as well. The median hardness ratio value for 
the CC clusters is ~ 0.85 while for NCC clusters it is ~ 1.0 in this 
range. 

the point sources removed, we then used the tool specex- 
tract to extract spectra with no grouping. The spectra 
include energies between 0.5 and 8 keV from the point 
source-cleaned level 2 event files in the defined regions for 
fitting in XSPEC. Since we need the bulk temperatures 
as a proof of principle and not for use with detailed clus- 
ter temperature analyses, if a cluster has more than one 
observation, we only utilized the longest exposure event 
files; using the spectrum from such a large region of each 
cluster provided us with more than enough counts for 
high S/N in each energy channel in the fitting process. 
Additionally, the spectra for front- and back-illuminated 
chips were extracted separately since they have different 
sensitivities. 

Appropriately subtracting the background spectrum is 
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vital if an accurate and believable cluster spectral fit is to 
be achieved. The same procedure was followed as for de- 
termining the background counts for the Chandra images 
discussed above, but we instead used the ACIS Blank- 
Sky Background files available from the CIAO database. 
Our temperature results below agree with published tem- 
peratures within our errors. All our spectral fits used 
the energy range (0.5 — 7) keV. The extracted event 
and normalized background spectra were then taken into 
XSPEC, where we used an APEC model with photoelec- 
tric absorption to fit the background subtracted cluster 
spectra. Weigh ted average values for ga l actic hydrogen 
colum n density (jDickev fe Lo ckman 1990; Kalberl a et alJ 
I2005T) were obtained from the online HEASARC Nh 
tool 9 but were allowed to change during fitting. Red- 
shifts were taken from Table [TJ and initial temperatures 
were from the M — T relation of Equation [3] 

The results for three paired clusters can be found in 
Table [6] The masses of the CC and NCC clusters in the 
first two matched pairs are within 10% of each other, 
while the third (unmatched) pair is given for further il- 
lustration. The regions over which the spectra were ex- 
tracted are provided, along with fit results for hydro- 
gen column, metallicity, and the fraction of the observed 
count rates attributable to the source instead of the back- 
ground. The reduced \ 2 an d degrees of freedom for each 
fit are also provided. Most fitted hydrogen column values 
are within a factor of two or three from those obtained by 
the HEASARC online tool, discrepancies attributable to 
the use of the ACIS blank-sky backgrounds for the back- 
ground subtractions. These backgrounds include data 
from several pointings and so likely include data from a 
region of the sky with slightly different galactic hydrogen 
than the regions of our particular cluster observations. 
A normalized temperature is also given, which is sim- 
ply the spectroscopically fit temper ature Tf a divided by 
the bulk cluster temperature Tchen (jChen et al.ll2007t ) in 
Table[TJ It is important to recall that Tchen has been cor- 
rected for the influence of cooler gas in CC cluster cores. 
In addition, recall that CIAO and CALDB 4.1.1 were 
used to update the analyses for the two mass-matched 
cluster pairs while the analyses for A4059 and A3391 
(slightly cooler clusters at 3.94 and 5.89 keV, respec- 
tively) were completed with CIAO and CALDB 3.3.0. 
Lastly, a fourth cluster pair was also identified, including 
the CC cluster A2204. Measured temperatures for A2204 
were exceedingly high, and a literature search revealed 
that these Chandra data over-estimate the temperature 
of A2204 compared t o the XMM and Suzaku satellites 
(|Reiprich et alJ I2008D , possibly caused by uncorrected 
PSF smearing in the XMM data and in accurate Chan- 
dra c alibration at high temperatures (see iReiprich et alJ 
120081 and references therein) . In light of this discovery, 
we removed A2204 and its pair from the bulk tempera- 
ture analysis. 

In general, an offset similar to that found in the ob- 
served normalized hardness ratios is present in the nor- 
malized temperatures. The observed normalized hard- 
ness ratio and temperature profiles are plotted for each 
paired cluster in Figure [TT1 In most cases, the tempera- 
ture profiles track the observed hardness profiles reason- 
ably well, indicating that the hardness ratio offset does 

9 http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 



indeed originate in an offset in cluster temperatures and 
lending further confidence to our assumption that hard- 
ness ratio maps can act as first-order proxies for clus- 
ter temperature. The temperature and hardness ratio 
profiles do not agree in every instance, but we believe 
the discrepancies are well explained by temperature sub- 
structure in some of the clusters (see below), which could 
easily be the result of minor cluster mergers. Relatively 
large and numerous pockets of temperature substructure 
create significant deviations in the average temperature 
profile. The hardness ratios do not manifest the substruc- 
ture as strongly owing to the hardness ratio's derivation 
from broad-band images, which tend to suppress temper- 
ature substructure due t o the n 2 dependence of surface 
brightness. Additionally, ICavagnolo et al.l (|2008l ) found 
that the ratio of hard to soft band spectroscopically fit- 
ted temperatures tends to be higher in systems known 
to be merging, where a deviation from unity indicates a 
band-dependence to the temperature fit. This is evidence 
suggesting unresolved clumps can alter the fitted temper- 
atures in a band-dependent fashion. In their study, the 
hard band ratios for A401 and A2255 were both con- 
siderably greater than 1.1 (1.37 and 1.32, respectively), 
which seems to corroborate our suggestion that tempera- 
ture substructure in our clusters explains the differences 
between their normalized hardness ratio and fitted tem- 
perature profiles. 

For each cluster, Figure [TTJ also contains ex- 
pected normalized hardness ratios corresponding to the 
spectroscopically-fit temperatures, which are plotted as 
dotted lines. These expected normalized hardness ra- 
tios were calculated using WcbSpec and the temperature 
results from fitting the extracted spectra with XSPEC. 
They include no metals since we are interested in prob- 
ing the effects of the temperature alone on the hardness 
ratios. The hardness ratios were normalized by the cor- 
responding hardness ratio for a cluster with zero metals 
exhibiting the same virial temperature as the cluster in 
question. As can be seen in Figure [Til m most cases 
the expected normalized hardness ratios agree well with 
both the observed normalized hardness ratios and the 
normalized spectral temperatures, except where signif- 
icant temperature substructure is apparent (e.g., A401 
see below). 

One should not be concerned with the increased un- 
evenness in the temperature profiles as compared to the 
relatively smooth hardness ratio profiles. Recall from 
Figure [7] that the hardness ratio has a clear but mod- 
est dependence on temperature. This has the effect of 
"flattening" the hardness ratio profile compared to the 
temperature profile from which it is derived, removing 
a great deal of the unevenness from pockets of cooler 
and/or warmer gas. Below, we discuss the remaining dis- 
crepancies between temperature and hardness ratio pro- 
files for each cluster separately with the aide of published 
temperature maps. 

4.3.1. A478 

The normalized temperature of A478 shows a large 
spike roughly three arc minutes from the cluster co re. A 
temperature map from lBourdin fc Mazzottal (|2008l ) con- 
tains a ring of hotter material between two and four ar- 
cminutes from the cluster core, ~ 2 — 3 keV warmer than 
the material just inside it. This is enough substructure to 
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explain the large spike in our temperature profile. But, 
no comparable substructure is seen in the X-ray surface 
brightness map. 

4.3.2. A401 

A temperature map from iBourdin fc Mazzottal (|2008f ) 
of A401 shows large blotches of gas both cooler and 
warmer than the average cluster temperature. There is 
one large warm region roughly two arcminutes in radius 
between three and six arcminutes from the core ~ 3 keV 
warmer than the average cluster temperature. Several 
cool pockets the same radial distance from the core as 
the warm blotch are present as well, but are only ~ 1.5 
keV cooler than the average temperature. It is possible 
radially averaging the temperature in this region results 
in an overall increase in temperature compared to the 
average, which is what is observed in the temperature 
profile for A401 in Figure QTJ relative to the hardness 
ratio profile. 

4.3.3. A85 

A ring of small warm gas clumps nearly twice the aver- 
ag e temperatur e of A85 is evident in a temperature map 
in lDurret et al.l (|200 5a). These small clumps appear be- 
tween three and seven arcminutes from the cluster core, 
and line up nicely with the slightly elevated tempera- 
tures in our temperature profile. Unfortunately, A85 is 
an example of the limitations of looking beyond the core 
region; we completely ignore the fairly relaxed central re- 
gion in hopes that cool core gas will not bias our results, 
yet we are confronted with more temperature substruc- 
ture. Again, there is no comparable substructure in the 
surface brightness and hardness ratio maps. 

4.3.4. A2255 

Our temperature profile for A2255 matches very well 
with the observed hardness ratio profile, although there is 
a small dip in temperature around t hree arcminutes from 
the co re. A temperature map from ISakelliou fc PonmanI 
(|2006| ) reveals a patch of cooler gas roughly three to four 
arcminutes from the core that accounts for the observed 
temperature dip. The temperature profile shows evi- 
dence it is beginning t o increase four or five arcm inutes 
from the core, and the ISakellio u fc Ponm au ( 200(3) tem- 
perature map does contain several significantly warmer 
patches at roughly this radial distance. 

4.3.5. A4059 

A4059 is another cluster whose normalized tempera- 
ture profile agrees well with the observed normalized 
hardness ratio profile. A slightly elevated temperature 
at roughly four arcminutes from the core is ob s erved , 
however. The temperature map of lRevnolds et alj (|2008D 
shows some patches of warmer gas a few arcminutes from 
the core, but it is hard to reliably correlate these warm 
gas pocke ts with the structure in our temperature profile 
since the iRevnolds et al.1 ()2008l ) temperature map does 
not have the same angular extent. Their map focuses on 
the inner cluster region, which we have excluded. 

4.3.6. A3391 

The only discrepancy between the hardness ratio 
and temperature profiles for A3391 is a rapid drop-off 



in temperature beyond six arcminutes. Although the 
resolution is poo r , an ASCA temperature map from 
iMarkevitch et all (|1998h reveals larger cooler regions 
beyond roughly six arcminutes, again matching our 
temperature results very well. 



Overall, the results of the hardness ratio and temper- 
ature analysis appear to indicate that there are both 
similarities and differences between CC and NCC clus- 
ters well outside the cluster cores. A statistically signifi- 
cant offset in normalized hardness ratios between CC and 
NCC clusters suggests that on average CC clusters tend 
to have more cool gas well beyond their cores than NCC 
clusters. In addition, there is general agreement between 
hardness ratio and spectral temperature profiles, further 
suggesting there is in fact a difference in ICM tempera- 
ture between CC and NCC clusters that drives the offset 
in hardness ratios. Although the hardness ratio profiles 
are quite fiat, cluster temperature maps reveal the pres- 
ence of temperature substructure, which is manifest in 
the much more irregular temperature profiles. Such tem- 
perature substructure is expected from our simulations in 
not only NCC but CC clusters as well since CC clusters 
suffer late minor mergers. In the next section, we com- 
pare these results to an analogous analysis of simulated 
clusters. 





Fig. 11. — Normalized hardness ratio and normalized spectro- 
scopically fit temperature profiles for the six mass-matched clus- 
ters. CC clusters are in blue while NCC clusters are in red, with 
mass-matched clusters in the same row. The thick solid lines with 
error bars are the observed hardness ratio values with 90% regions. 
The two solid lines delineate the 90% regions for the temperatures, 
which have been normalized by the cluster temperatures in Table 
[T] Dotted lines correspond to expected model hardness ratios for 
clusters with 0.3 solar metallicity and the redshifts and temper- 
atures from Table [6] The first four clusters were analyzed with 
CIAO and CALDB 4.1.1, while the last two clusters were analyzed 
with CIAO 3.3.0. 
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5. SIMULATION COMPARISONS 

With the analysis of our observed cluster sample in 
hand, we next compare these results to those of the sim- 
ulated clusters from B08. What can we say about the 
excess cool gas found in CC clusters outside their cores? 
What are the similarities and differences between the ob- 
served and simulated properties of CC and NCC clusters? 
How might we improve our simulations so that they bet- 
ter represent the physics of observed clusters? We com- 
pare the results of the simulated clusters to those of the 
real clusters to answer these questions and learn how to 
achieve more realistic clusters in future simulations. We 
begin by providing a quick description of the simulations 
from B08, then discuss the results of the surface bright- 
ness profile analysis, and conclude with an overview of 
simulated temperature and hardness ratio profiles. 



5.1. Simulation Description 

The simulations w ere performed using the Enzo code 
(jO'Shea et all I2005D . This code combines an N-body 
algorithm that evolves collisionless dark matter parti- 
cles using an Eulerian hydrodynamics scheme (Picccwise 
Parabolic Method or PP M) that employes adaptive mesh 
refinement (AMR) (see IBurns et all 120041: iMotl et all 
l2004t IBurns et all [20081 ). These simulations include ra- 
diative cooling, star formation and energy feedback. Ra- 
diative cooling i s based on a Raymon d-Smith plasma 
emission model (jBrickhouse et al1ll995() that assumes a 
constant metallicity of 0.3 solar. Star formation helps 
soften the cool cores by allowing excessively cool gas to 
drop out into "star particles." These star particles also 
allow energy to be reinserted into the sys tem by means of 
supern ovae by a prescrip tion outlined in lCen fc Ostrike"rl 
(|1992f) and IBurns et all (|2004f ). This is explained in de- 
tail in B08. 

In B08, we examined the simulation volumes and suc- 
cessfully compared bulk X-ray properties to those ob- 
served in surveys of galaxy clusters. For the first time, 
our simulations produced both CC and NCC clusters in 
the same numerical volume. B08 used comparisons of 
spectroscopic-like temperatures, f gas (the fraction of gas 
mass over total mass), M^oo,Gas (the gas mass within 
rsoo) and indi vidual /3-mq d el fits to similar observational 
collections in I Chen et al.l (|2007t ) . We showed that our 
simulations match observed bulk properties of clusters 
well enough for use in making testable predictions that 
we examine in this paper. 

CC and NCC clusters experience different merger his- 
tories according to our simulations. NCC clusters suffer 
early major mergers that destroy any cool cores that have 
begun to form, while CC clusters remain undisturbed 
during early epochs. It is important to note that star 
formation and feedback play a crucial role in "soften- 
ing" the cool cores (i.e., heating, expanding, and lowering 
the gas density), thus allowing collisions with compara- 
ble mass halos to ram pressure strip and destroy nascent 
cool cores. At later times, CC clusters may be involved 
in mergers, but their cool cores are large and "strong" 
enough to not be disrupted by the event. In addition, 
fewer larger mergers in the histories of CC clusters means 
that the gas through the clusters is shock-heated less and 
so remains cooler than gas in NCC clusters. 



5.2. Simulated Surface Brightness Profiles 

The azimuthally averaged surface brightness profiles 
of 10 CC and 78 NCC simulated clusters, each normal- 
ized to its expected luminosity (M 3 ' 2 ), are displayed in 
Figure [T2l with 90% regions. This is a statistically com- 
plete sample that includes every cluster at z = with 
M > 10 14 M that can be clearly identified as a CC or 
NCC cluster and that lack extreme recent merger effects; 
see B08 for more details. Double and single /3-models 
were fitted to the averaged CC and NCC cluster profiles, 
respectively, and over-plotted. The /3-model fits to the 
simulations indicate that fitting observed surface bright- 
ness profiles to the farthest edge of observations (for 
Chandra, often to ~ 0.3 — 0.5r2oo) may overestimate the 
surface brightness at large radii, but more so in CC than 
NCC clusters. This overcstimation of surface brightness 
in turn leads to an overestimation of th e extrapolated 
cluste r gas mass at r2oo as discussed in Hallman et al.l 
(|2006h . Clearly, effects such as this complicate the use of 
galaxy clusters for precision cosmology and suggest that 
care must be taken when experiments are based on their 
extrapolated properties. 

Figure [T3] shows the averaged observational data com- 
pared to the averaged surface brightness profiles of our 
simulated clusters; 90% regions have been removed for 
clarity. While the simulated CC clusters show an ex- 
cess luminosity at the core relative to observations, the 
overall agreement between simulations and observations 
is reasonable, though it is clear more work is required 
to better model observed clusters. The CC core excess 
luminosity is a known limitation due to the overcooling 
of high density regions within the simulations and is dis- 
cussed in B08. Also, the larger errors associated with 
the simulations at small radii in Figure [T2l are a result of 
limited resolution for the simulated clusters (16 kpc peak 
resolution). Pixel sizes at the simulated cluster centers 
are relatively large and provide only a few points with 
which to average across, resulting in 90% regions that 
appear to underestimate the precision of the simulated 
data below ~ 0.1 r2oo- 

As was done for the observations, double and single (3- 
modcls for CC and NCC clusters, respectively, were fitted 
to the cluster profiles individually between and 0.4r2oo- 
These individual fits were then grouped into cluster cat- 
egory, averaged and displayed in Table Histograms 
of the individual cluster fits are provided in Figure 1141 
Much can be said by comparing the observed and simu- 
lated results in Tableland Figurc[T31 and there are some 
clear similarities. K-S test results confirm that there is no 
statistical differences between outer CC and NCC r corey 2 
fits, neither for observed nor simulated clusters, a result 
that co ntrasts with those found using solely si ngle (3- 
models (IVikhlinin et al.l l2006t IBurns et al.l I2008D where 
the cool core is excluded from the fit. Additionally, the 
size of the averaged observed and simulated CC cluster 
cores (r corej i) are consistent. The profiles also clearly 
indicate that simulated NCC clusters do a good job of 
capturing the flat inner profile inside cluster cores. 

There are important differences in /3-model parame- 
ters and surface brightness profiles, however. The effect 
of overcooling the simulated cluster cores is to drastically 
raise the inner /3 slope for the CC clusters compared to 
the observed CC clusters. In the outer profiles, and as 
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evident in the (3 slopes in Table [2 the observed clusters 
exhibit flatter profiles than the simulations. This may 
be evidence for additional heating and feedback not cap- 
tured by our "star particle" technique in the outer regions 
of observed clusters, perhaps puffing up the outer profile 
and making them brighter - another clue that our feed- 
back prescriptions must be updated. This also suggests 
that the entropy in the outer regions of observed clusters 
may be higher than in our simulated clusters. All of these 
differences point to corrections that must be made in fu- 
ture iterations of the simulations if we wish to accurately 
model observed clusters. 
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Fig. 14. — Histograms of simulated /3-model fit parameter val- 
ues. (Top) Outer r core fits. (Bottom) /3 fits. It is clear from the 
histograms that both r CO re and j3 fits for CC (blue) and NCC (red) 
clusters are statistically indistinguishable. 
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Fig. 12. — The averaged surface brightness profiles of 10 CC 
(blue) and 78 NCC (red) simulated clusters. The data have been 
averaged across all the clusters of each type and uncertainties are 
displayed as 90% regions demarcated as dotted lines. A /3-model 
has been fit to the region — 0.4r200 separately for both cluster 
types and then extrapolated over the full range of the plot. The 
vertical dashed line indicates the edge of the fitting region. The 
NCC profile has been arbitrarily rescaled by a factor of 10 to sep- 
arate the curves for clarity. The /3-model somewhat overestimates 
the surface brightness at the virial radius, particularly in CC clus- 
ters. 
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Fig. 13. — A comparison of simulated and observed averaged 
surface brightness profiles. The simulated clusters include 10 CC 
(blue) and 78 NCC (red), all the available clusters at z = 0. The 
NCC clusters have been arbitrarily divided by 10 to separate the 
lines for clarity. Uncertainties have been removed for clarity. The 
excess central luminosity of the simulated CC clusters is due to 
overcooling as discussed in B08. The outer profiles for the simu- 
lated clusters, especially the cool cores, are somewhat steeper than 
the observed clusters (see Table [2}. 



5.3. Simulated Temperatures and Hardness Ratios 



Unlike for the observed clusters, we have full knowl- 
edge of simulated temperatures and so we calculated av- 
eraged temperature profiles for the simulated clusters. 
These give further confidence in the use of observed hard- 
ness ratios as a first-order proxy for cluster temperatures, 
and provide another test of the assumption that the off- 
set seen between the observed and simulated CC and 
NCC hardness ratio profiles is a consequence of differing 
cluster temperatures beyond the cores of mass-matched 
cluster pairs. Figure [15] shows the averaged tempera- 
ture profiles, again with standard deviation error bars. 
It is clear in the figure that the simulated temperature 
profiles exhibit the same offset seen in the hardness ra- 
tio profiles one would expect if in fact the temperatures 
drive the hardness ratio values; CC clusters are cooler 
and NCC clusters are warmer well beyond the cluster 
cores. Another K-S test over the region 0.3 — 1.0r2oo for 
the averaged temperature profiles yields a probability of 
8.4% that the temperature profiles come from the same 
distribution. This result is compelling evidence that the 
offset does indeed exist in temperature. 

We also performed a hardness ratio analysis on the 
simulated clusters for comparison with observations. Fol- 
lowing the same procedure for simulated data as we used 
with the Chandra data (detailed in Section [4. 2p . we split 
the synthetic X-ray images of the simulated clusters into 
two energy bands, 2-8 keV and 0.5-2 keV. We then di- 
vided the two images (2-8 keV/0.5-2 keV) and normal- 
ized them to the hardness ratio that corresponds to the 
"virial" temperature to produce normalized hardness ra- 
tio images. As in the observed hardness ratio maps, pixel 
values less than one represent gas with T < T V i r and pixel 
values greater than one represent T > T v i r . Also, as in 
the observed clusters an M—T relation is used to find the 
cluster temperatures for normalization. Plotting the log 
of simulated cluster masses M500 versus log of the clus- 
ter temperatures, a best-fit line yields the corresponding 
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mass-temperature relation: 

l0g - (5^^)-°- 629 + 1 - 551 ° g -(4i^v)- 

(5) 

The observed relation from IChen et al.l (|2007f ) and the 
simulated cluster relation agree well, where both are con- 
sistent with the 3/2 slope expected if the clusters arc 
scale-free. In an effort to mimic what is done with ob- 
servations as much as possible, we "corrected" the simu- 
lated CC cluster temperatures for the presence of their 
cool cores. We first fit for the M — T relation for sim- 
ulated NCC clusters only and found the deviation from 
this relation for each CC cluster separately We then fit a 
line to the resulting temperature offset-cluster mass plot 
to define an average temperature correction for a given 
mass. Finally, we applied this temperature correction to 
all simulated CC clusters. 

As we did for the observed clusters, simulated clusters 
were mass-matched to within 10%. The hardness ratio 
profiles for the individual simulated clusters are found 
in the top of Figure HH Mass-matched pairs are plotted 
with the same line style and line thickness. The bot- 
tom of Figure [16] contains the corresponding averaged 
CC and NCC hardness ratio profiles. In the outer re- 
gions of the clusters, we qualitatively see the same hard- 
ness ratio offset as was observed in the real cluster pro- 
files: the NCC clusters are harder than the CC clusters, 
despite very similar masses. Performing the same K-S 
test as was done for the observed clusters, but for the 
region (0.3 — 1.0)?'2oo to avoid the region of simulated 
clusters most likely affected by the overcooled cores, the 
outcome was identical. Namely, the probability that the 
simulated CC and NCC hardness ratio pixel sets come 
from the same distribution is <gc 1%. However, there 
are also differences in the hardness ratio profiles between 
observations and simulations. In comparing Figure [5] to 
Figure [TBI the observed clusters have a larger separation 
in hardness ratio between CC and NCC clusters in the 
radii 0.1 to 0.4 r2oo- 
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Fig. 15. — Averaged simulated temperature profiles for CC (blue) 
and NCC (red) clusters, normalized by the M — T relation tem- 
peratures. To mimic what is done for observations, the CC cluster 
temperatures have been "corrected" for their cool cores. Error bars 
delineate the standard deviation at each radial point. 

6. SUMMARY AND CONCLUSIONS 

Using an optimized data pipeline to reduce Chandra 
and ROSAT data, we analyzed 30 nearby, rich Abell 




Fig. 16. — (Top) Individual simulated hardness ratio radial pro- 
files for CC (blue) and NCC (red) clusters. The mass-matched 
pairs are next to each other in the legend and have the same line 
style and line thickness. (Bottom) The average hardness ratio pro- 
files. Error bars are standard deviations at each point. 

galaxy clusters split nearly equally between CC and NCC 
clusters. While the expected increase in radial informa- 
tion from the ROSAT data rarely appeared, the addi- 
tional data did help reaffirm the data reduction pipeline 
by matching well across instruments. The observational 
data were analyzed by examining radial surface bright- 
ness profiles and hardness ratio maps. 

Utilizing the ENZO AMR simulations, we also ana- 
lyzed 88 simulated clusters at z = 0. To the best of 
our knowledge, these are the first simulations to produce 
both CC and NCC clusters within the same cosmologi- 
cal volume. These simulations were used to create data 
products that are directly comparable to observed clus- 
ters and initial efforts indicate that the simulations are 
reasonable first-order representations of actual clusters, 
but several modifications must be made for the future. 
The simulations were then used to make testable predic- 
tions about the differences between CC and NCC clus- 
ters. Specifically, simulations suggest that NCC clusters 
suffer early major mergers (roughly defined as > 50% in- 
crease in mass in less than one Gyr) that disrupted any 
nascent cool core that may have begun to form. CC clus- 
ters, on the other hand, settle down and strengthen for 
longer periods of time and suffer generally smaller merg- 
ers; this allows their cores to survive subsequent accretion 
events. 

Our major results include 

• Observed CC and NCC clusters, fit with single 
and double /3-models, respectively, have compara- 
ble cluster core radii, in contrast to that claimed 
from single /3-model fits where cool cores were ex- 
cluded. Apparently, this technique for CC clusters 
is not completely successful and biases the fits to 
smaller core radii. 

• There is a statistically significant difference in the 
slopes (i.e. j3 values) of the outer surface brightness 
profiles between observed CC and NCC clusters. 
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This may be consistent with a different evolution- 
ary history between these cluster types. 

• CC clusters show a significant offset in normalized 
hardness ratio values beyond their cool cores in 
comparison to NCC clusters, which is consistent 
with cooler gas extending out to at least 0.4r2oo hi 
CC clusters. 

• Simulated clusters have steeper outer surface 
brightness profile slopes (/3 values) than observed 
clusters (Table El). This may suggest some missing- 
physics in the simulations beyond the cluster core 
(i.e., preheating?). 

Although our simulations generally agree qualitatively 
with observations, comparisons with observations also 
suggest a number of inconsistencies that must be resolved 
if we desire to use simulated clusters as statistical proxies 
for real clusters in precision cosmological studies. Inade- 
quacies include overcooling in the CC cores, overly simi- 
lar and overly large (3 slopes between simulated CC and 
NCC clusters at intermediate cluster radii (~ 0.4r2oo) 
compared to observed clusters, and a smaller hardness 
ratio offset between CC and NCC clusters than found in 
observations. There is also evidence that observed clus- 
ters may have higher entropy than our simulated clusters 
via their fla tter outer su r face b rightness profiles. 

Recently, lArieli et"aT1 (|2008f ) suggest the use of Gal- 
cons, or galaxy constructs, in cosmological simulations 
instead of " s tar p articles" for regulating heat feedback. 
lArieli et alj (|2008f ) find that star particles tend to group 
in dense and cold regions of the ICM gas, namely, the 
cores of CC clusters where the gas conditions arc more 
favorable to star formation. Any feedback heating from 
these star particles has little effect on the cluster beyond 
the inner regions of the core where the stars form. Gal- 
cons, however, are capable of feedback heating over much 
more of the cluster in the form of "galactic winds." These 
winds can fill the role of any number of proposed feed- 
back/heating mechanisms that may be at work in real 
clusters. The use of Galcons in our simulations has the 
potential to ameliorate our CC cluster overcooling prob- 
lem while at the same time increasing the effectiveness of 
feedback heating in in-falling halos, which could flatten 
the overly steep surface brightness profiles, and perhaps 
other unexpected yet positive consequences. 

We plan to extend this direction of research with addi- 
tional simulations that include Galcons and AGNs (e.g. 
iHoeft et all 120081 iPuchwein et~ai1 l2008f ) to take advan- 
tage of potential AGN feedback mechanisms. In addi- 
tion, the clusters selected for this study were chosen be- 
cause they all have at least 3 kscc of observations with 
XMM, so we plan to include XMM data in our obser- 
vational analysis as well. With additional observations 
and new simulations, we hope to further understand the 
differences between CC and NCC clusters and their use 
in precision cosmological studies. 
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TABLE 1 

Observed Cluster. Properties and Pointing Data 



Cluster 


(C)handra Obs ID 
(R)OSAT Obs ID 


ACIS-I/S a 
PSPC-A/B a 


Exp 
Time (ksec) 




Tchen 

(keV) 


M 2 oo c > d 
(10 14 Ms) 


r200 c 
(Mpc) 


Cool b 
Core 


Previous 
Publications 


A85 


(C) 904 
(R) 174A 
(ft) 174A1 

111/ -LI -L 

(R) 250 


I 

B 
]j 

B 


39 
2 
3 

10 


0.0521 


6.51 


8.27 


1.71 


Yes 


1, 2 


A119 


(C) 4108 
(R) 251 


9 
B 


10 
15 


0.0444 


5.69 


9.12 


1.67 


No 


3, 4 


A262 


(C) 2215 
(R) 254 


S 
B 


29 
9 


0.0163 


2.25 


0.93 


1.00 


Yes 


5, 6, 7 


A401 

A (OCl 1 


(C) 518 
(C) 2309 
(R) 182 
(R) 235 


I 
I 

B 
B 
i 


18 
12 
7 
7 


0.0748 


7.19 


8.77 


2.01 


No 


8, 9, 10 



A1656 



A1795 



(C) 503 
(C) 1513 
(C) 3209 
(C) 4289 
(C) 4946 
(C) 4947 
(R) 32A 
(R) 32A1 

(R) 33 
(R) 34A 
(R) 34A1 
(R) 35A 
(R) 35A1 
(R) 186 
(R) 521 



(C) 555 
(C) 556 
(C) 1086 
(C) 1112 
(C) 1113 
(C) 1114 
(R) 5 
(R) 6 
(R) 9 
(R) 13 
(C) 493 
(C) 494 
(C) 5289 
(R) 55 
(R) 105 
(R) 145 
[Cj 891 
(C) 4977 
(C) 6101 
(R) 161 
(R) 249 



S 
S 

s 
s 
s 
s 

B 
B 
B 
B 
B 
B 
B 
B 
B 



I 

S 

s 
I 
I 
I 

B 
B 
B 
B 

IT 
S 
I 

B 
B 
B 

IT 
S 
I 

B 
B 



9 
25 
97 
97 
21 
30 

1 

31 
15 
10 
21 

6 

26 
5 
3 



A478 


(C) 1669 
(C) 6102 
(C) 7217 
(R) 193 


S 
I 
1 

B 


42 
10 
19 
22 


0.0882 


6.91 


9.38 


1.82 


Yes 


13, 14 


A496 


(C) 931 
(C) 3361 
(C) 4976 

(R) 24 


S 

s 
s 

B 


19 
10 
76 

9 


0.0331 


3.91 


4.83 


1.35 


Yes 


15, 16 


A1060 


(C) 2220 
(R) 135 
(R) 200 


I 

B 
B 


32 
20 
16 


0.0124 


3.15 


2.46 


1.20 


Yes 


17, 18 


A1367 


(C) 514 
(R) 153 


S 
B 


41 
19 


0.0214 


3.55 


7.36 


1.29 


No 


19, 20 



9 
10 
10 
10 
10 
10 
21 
22 
20 
21 



0.0231 8.07 



13.7 



2.00 No 21, 22 



A2029 



20 
20 
15 
26 
36 
18 



0.0622 6.17 



10.21 



1.60 



6, 23 



A2052 



(C) 890 
(C) 5807 
(R) 275 



20 
79 
10 

3 

13 



0.0766 7.93 



10.43 



2.08 



24, 25 



S 
S 
B 



37 
129 

6 



0.0353 3.12 



2.72 



1.18 



Yes 



26, 27 



A2063 


(C) 4187 


I 


9 


0.0355 


3.56 


2.38 


1.31 


Yes 


28, 29 




(C) 5795 


S 


10 
















(C) 6262 


S 


14 
















(C) 6263 


s 


17 
















(R) 128 


B 


10 
















(R) 184A1 


B 


10 
















(R) 376 


B 


7 














A2142 


(C) 1196 


S 


12 


0.0894 


8.46 


15.21 


2.06 


Yes 


30, 31 




(C) 1228 


S 


12 
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TABLE 1 — Continued 



Cluster 


(C)handra Obs ID 
(R)OSAT Obs ID 


ACIS-I/S a 
PSPC-A/B a 


Exp 
Time (ksec) 




Tchen 

(keV) 


M 200 c ' d 
(10 14 Mq) 


r200 c 
(Mpc) 


Cool b 
Core 


Previous 
Publications 




(C) 5005 
(R) 96 
(R) 233 
(R) 415 
(R) 551 


I 

B 
li 
13 
B 


45 
6 
5 

19 
6 














A2199 


(U) 497 
(C) 498 
(R) 644 


9 

S 

B 


20 
19 
41 


0.0299 


4.28 


4.29 


1.40 


Yes 


32, 33 


A2204 


(C) 499 
(C) 6104 
(R) 281 


9 
1 

B 


10 
10 

5 


0.1524 


6.38 


6.53 


1.87 


Yes 


34, 35 


A2255 


(C) 894 
(R) 512 


I 

B 


40 
15 


0.0809 


5.92 


8.28 


1.82 


No 


36, 37 


A2256 


(C) 965 
(C) 1386 
(C) 1521 
(C) 2419 
(R) 162A 
(R) 162A1 
(R) 163 
(R) 339 
(R) 340 
(R) 341 


9 
I 

s 
s 

B 
B 
B 
B 
B 
B 


11 

13 

3 
12 

1 

5 

11 

5 
9 

10 


0.0581 


6.83 


12.48 


1.91 


No 


38, 39 


A2319* 


(C) 3231 

(R) 73 
(R) 73A1 


I 

B 
B 


15 
1 
3 


0.0555 


9.12 


18.6 


2.12 


No 


40, 41 


A2597 


(U) 922 
(C) 6934 
(C) 7329 
(R) 112 


9 

s 
s 

B 


40 
53 
61 
7 


0.0852 


4.2 


3.92 


1.30 


Yes 


42, 43 


A3158 


(C) 3201 
(C) 3712 
(R) 310 


I 
1 

B 


25 
31 

3 


0.0590 


5.41 


5.93 


1.68 


No 


33 


A3266 


(C) 899 
(R) 211 
(R) 552 


I 

B 
B 


30 
7 

14 


0.0545 


7.72 


19.75 


1.94 


No 


44, 45 


A3391 


(C) 4943 
(R) 79 
(R) 80 


I 

B 
B 


19 
3 
7 


0.0550 


5.89 


6.20 


1.65 


No 


29 


A3558 


(C) 1646 
(R) 76 


S 
B 


15 
30 


0.0477 


5.37 


6.84 


1.66 


No 


29, 46 


A3562 


(C) 899 
(R) 237 


I 
B 


30 
20 


0.0502 


5.16 


3.59 


1.57 


No 


47, 48 


A3571 


(C) 4203 
(R) 287 


S 
B 


34 
6 


0.0397 


6.80 


8.86 


1.83 


Yes 


49, 50 


A3667 


(C) 5751 
(C) 5752 
(C) 5753 
(R) 234 


I 
I 
1 

B 


130 
61 
105 
13 


0.0530 


6.28 


5.41 


1.85 


No 


51 


A4038 


(C) 4188 
(C) 4992 
(R) 354 


I 
1 

B 


6 

31 

3 


0.0281 


3.22 


2.58 


1.21 


Yes 


52 


A4059 


(C) 5785 
(R) 175 


S 
B 


93 
5 


0.0456 


3.94 


4.49 


1.39 


Yes 


53, 54 
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17 ISato ct al. (2007), 18 Havaka wa et all 



2000 


), 28: |Kanov et al. 


120061)" 


2006 


), 33|L'o"*"*"r*e"t**"n 1 


20061). 



12009) 
120071) 



id 1200*51) 

24 IClarke et al.l l2005al ). 25lMivoshi et al.l T2005I). 26lBlanton et al 



120071), 35ISanders et al.l 



34^^^Hrfiet^d 

38IUlarke Hz tlinsslinl 1120061). 3 !J IMiller et al.l 12003*1. 40lO'Hara et al.l J20Q. 
12005b[), 44IFinoguenov et al.l 120061), 45ISauvageot et al.l 120051), 46[P*****setti et al 
49 lHudaverdi et al.l 120051) , 50 rVenturi et al.l 120021) , 51 IBriel etall 120041) , 52 lSlee 



2 9IDunn ^TT abian (200(g), 30ITittlev fc"* Henrikscn 1200^.,, 

et al.l 120051), 36lGovoni et al.l 12006*) 
304T). 41IGovoni et al.l 120041). 42lMorris 

' ' 120071), 47 : Finoguciio v et al.l 12004D, 4 8 TVenturi et al l 120031 ), 
Rovl 119981) . 53 lRevnolds et al.l l2008) , 54 C hoi et al.f 120041) 



12007 ), 27 Rizz a et al.l 
31IMarkevitch et al.l 120001), 32ISanders fc Fabiaa 
37ISakelliou fc Ponmanl 120061), 
"Fabianl 120051). 43lClarke et al.l 



Note. — Observed X-ray properties and instrument pointing data for the selected Abcll clusters. 

Chandra or ROSAT detector used. [p'Hara ct al. (2006]) c Calculated from scaling relations, sec text. O'Hara ct al. (2006), Finogucnov ct al. 



120011) , [Navarro et all U997D , IVoitl l2005l) d 
from IChen et aT] 120071) . 



Afsoo and Tchen taken from Chen ct al. (2007) and Ikcbc ct al. (2002)° Mass and temperatures not 
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TABLE 2 
Average /3-Model Fit Results 





Average f3\ 


Average r corej i 


Average $2 


Average r corEi 2 


Observed CC clusters 
Observed NCC clusters 
K-S Test Probabilities 


0.54±0.07 


0.01±0.004 


0.55±0.04 
0.49±0.02 
0.002 


0.08±0.02 
0.08±0.01 
0.209 


Simulated CC clusters 
Simulated NCC clusters 
K-S Test Probabilities 


0.98±0.41 


O.OlrhO.Ol 


0.70±0.10 
0.71±0.22 
0.650 


O.lOrbO.01 
0.13±0.06 
0.277 



Note. — The averaged values of individual /3-modcls fit between and 0.4r2oo f° r observed and simulated clusters, along with the standard 
deviation of each population. (Sec 35.21 for more details on the simulated profiles). Subscripts 1 and 2 denote the averages of the innermost and 
outermost r core and (3 values, respectively. The r core fits arc in units of r200- K-S probabilities arc probabilities that the fits for CC and NCC 
clusters come from the same parent distribution. 

TABLE 3 

Observed Cluster /3-Model Fit Parameters 



ClustrGr 


!■■' 1 


' core,l 




'core, 2 


DOF 


Reduced 


A262 


0.738 


0.003 


0.503 


0.051 


14 


0.082 


A426 


0.396 


0.001 


23.468 


0.670 


14 


0.148 


A478 


0.406 


0.006 


14.937 


0.573 


11 


0.149 


A496 


0.478 


0.007 


0.489 


0.054 


14 


0.137 


A1060 


0.494 


0.008 


0.463 


0.049 


14 


0.130 


A1795 


0.488 


0.009 


0.487 


0.058 


11 


0.139 


A2029 


0.504 


0.012 


0.515 


0.070 


14 


0.138 


A2052 


0.514 


0.013 


0.542 


0.084 


14 


0.148 


A2063 


0.513 


0.014 


0.557 


0.088 


14 


0.135 


A2142 


0.512 


0.014 


0.578 


0.094 


11 


0.144 


A2199 


0.520 


0.014 


0.579 


0.090 


11 


0.130 


A2204 


0.523 


0.013 


0.575 


0.086 


14 


0.139 


A2597 


0.531 


0.014 


0.575 


0.087 


11 


0.135 


A3571 


0.539 


0.015 


0.594 


0.093 


11 


0.126 


A4038 


0.536 


0.016 


0.589 


0.089 


14 


0.107 


A4059 


0.537 


0.016 


0.601 


0.092 


14 


0.109 


A401 






0.498 


0.082 


17 


0.065 


A1367 






0.448 


0.078 


17 


0.065 


A1656 






0.488 


0.074 


17 


0.062 


A2255 






0.472 


0.073 


17 


0.066 


A2256 






0.499 


0.084 


17 


0.052 


A3158 






0.520 


0.090 


17 


0.057 


A2319 






0.503 


0.082 


17 


0.052 


A3266 






0.505 


0.082 


17 


0.054 


A3391 






0.495 


0.081 


17 


0.064 


A3558 






0.491 


0.078 


17 


0.144 


A3562 






0.487 


0.074 


17 


0.205 


A3667 






0.487 


0.075 


17 


0.197 



Note. — Individual fit values and statistics for first the observed CC and then NCC clusters. "1" represents the inner 
/3-model fits (for CC clusters only) while "2" represents the outer fits. DOF show the degrees of freedom in the fits. The fits 
for A426 and A478 are clear outliers and so are not included in the average values found in Table [2] 



TABLE 4 

Observed Cluster Spectral Fitting Results 



Cluster CC/NCC Region (r 2 oo) T Che „ a T fit b T) it /T Chen c N H (10 22 cm- 2 ) d z e Z(Z Q Y Reduced x 2 DOFs Source % h 

A478 CC 0.14- 0.18 6.91 7.63±o;;^ 1.10 0.277 0.0882 0.211 1.011 52 67.5 

0.18 - 0.22 6.91 7.04±q'.42 L02 °- 279 0.0882 0.189 1.356 52 57.4 

0.22 - 0.26 6.91 5.9lto;|@ 0.85 0.274 0.0882 0.212 1.530 52 47.6 

0.26 - 0.30 6.91 5.09to 8° 0.74 0.308 0.0882 0.347 0.558 52 38.7 

A401 NCC 0.11 - 0.1575 7.19 6.95±o;op 0.97 0.113 0.0748 0.355 2.017 52 90.5 

0.1575 - 0.205 7.19 6. 53to.H 0.91 0.123 0.0748 0.353 1.671 52 85.7 

0.205 - 0.2525 7.19 6.65±°;|g 0.92 0.084 0.0748 0.353 1.237 52 80.5 

0.2525 - 0.30 7.19 6.53±g;s| 0.91 0.121 0.0748 0.306 1.247 52 74.7 

A85 CC 0.12 - 0.165 (151 Km^^ 097 0.008 0.0521 0.351 0.967 52 88 

0.165 - 0.21 6.51 6.32±g;42 °- 97 °- 013 0.0521 0.355 1.574 52 81.3 

0.21 - 0.255 6.51 6.20tg 43 0.95 0.033 0.0521 0.148 1.341 52 73.4 

0.255 - 0.30 6.51 5.86±g % 0.90 0.010 0.0521 0.240 0.997 52 69.9 

A2255 NCC 0.11 - 0.1575 5.92 5.641°;^ 0.95 0.025 0.0809 0.300 1.208 52 83 

0.1575 - 0.205 5.92 ^■^^fi 0.92 0.028 0.0809 0.402 1.097 52 79.3 

0.205 - 0.2525 5.92 5.34±^;|g 0.90 0.021 0.0809 0.342 1.331 52 74.9 

0.2525 - 0.30 5.92 lm °- 000 0.0809 0.269 1.400 52 70.5 

A4059 CC 0.1 - 0.15 3S)4 4.12+%H L05 0.028 0.0456 0.605 1.019 52 87~3 

0.15 - 0.20 3.94 4.28±qH 1.09 0.021 0.0456 0.381 1.334 52 81.3 |jj 

0.20 - 0.25 3.94 3.86lg;3 4 0.98 0.055 0.0456 0.142 1.457 52 68.3 g 

0.25 - 0.30 3.94 4.26to.lt lm °- 061 0.0456 0.299 1.827 52 53.8 |. 

A3391 NCC 0.11 - 0.1575 5.89 5.47±q 7 5 \ 0.93 0.059 0.0550 0.043 0.968 52 85.9 

0.1575 - 0.205 5.89 5.18±g-J| 0.88 0.039 0.0550 0.000 1.070 52 81.1 S- 

0.205 - 0.2525 5.89 5.68lo;g7 0.96 0.046 0.0550 0.154 1.370 52 73.8 PL 

0.2525 - 0.30 5.89 ^Oto.lt °- 70 °- 104 0.0550 0.063 1.259 52 65.6 



NOTE. — Mass-matched CC/NCC pairs are separated by horizontal lines 



been corrected for cool cores. Fit results for galactic hydrogen column density. Rcdshifts from RTrlara ct aT| Fit results for cluster mctallicity. Degrees of freedom in the fit. Percent of extracted spectra counts 



attributable to the cluster in the corresponding spectral region 



